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Abstract. In a two-mode approximation, Bose-Einstein condensates (BEC) in a 
double-well potential can be described by a many particle Hamiltonian of Bose- 
Hubbard type. We focus on such a BEC whose interatomic interaction strength is 
modulated periodically by <5-kicks which represents a realization of a kicked top. In 
the (classical) mean-field approximation it provides a rich mixed phase space dynamics 
with regular and chaotic regions. By increasing the kick-strength a bifurcation leads to 
the appearance of self-trapping states localized on regular islands. This self-trapping 
is also found for the many particle system, however in general suppressed by coherent 
many particle tunneling oscillations. The tunneling time can be calculated from 
the quasi-energy splitting of the corresponding Floquet states. By varying the kick- 
strength these quasi-energy levels undergo both avoided and even actual crossings. 
Therefore stimulation or complete destruction of tunneling can be observed for this 
many particle system. 
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1. Introduction 

After their first experimental realization [1,2], Bose-Einstein condensates (BEC) have 
stimulated an enormous amount of theoretical investigations. Especially the recent 
progress in confining and manipulating BECs [3,4] affords an opportunity to realize a 
variety of quantum mechanical models experimentally 

Since a full many particle treatment of BECs is only possible for a small number of 
atoms, such systems are mostly described in the celebrated mean-field approximation, 
which describes the system quite well for large particle numbers at low temperatures. 
Most previous studies focused on relatively simple models like the two-mode system 
in order to investigate the correspondence between many particle and mean-field 
description [5-7]. Most interestingly, the mean-field description as a large particle 
number limit of many particle systems is formally related to the usual classical limit 
of quantum mechanics. In a number of recent papers consequences of the classical 
nature of the mean-field approximation are discussed and semiclassical aspects are 
introduced [8-12]. For a two-mode system even the eigenenergies and eigenstates of the 
many particle system could be reconstructed approximately from the mean-field system 
in a WKB type manner with astonishing accuracy [12]. Recent experiments with a 
relatively small number of particles [3] offer the opportunity to study the interplay of 
classical and genuine quantum behavior. 

Within the theory of quantum chaos systems with a periodic time dependence play 
an important role. Despite their apparently academic character systems with a delta- 
type time dependence are of particular relevance. Their prominence is not only due 
to their relatively simple handling in theory, but first of all to their typical behavior. 
Therefore systems like the kicked rotor or the kicked top became standard models in 
this area [13, 14]. 

In the present paper we focus on an iV-particle BEC in a double-well trap with a 
periodically kicked interaction strength. Experimentally the particle interaction may be 
modulated via a Feshbach resonance, a simpler realization, however, can be achieved by 
kicking the single particle tunneling coupling instead by varying the optical potential. 
As we will see later, these two systems are equivalent in the case of a symmetric double- 
well. In a two-mode approximation such a system can be described by a two-site Bose- 
Hubbard type Hamiltonian. In the case of a symmetric trap this system is in fact a 
realization of a kicked top studied, e.g., by Haake et al. in the context of quantum 
chaos [14,15]. Considering this BEC context, the question if the kicked top could be 
realized experimentally [16] can be answered in a novel way which offers the opportunity 
to study aspects of the system which were not investigated in the past. 

One of the most prominent features of the corresponding time independent system 
is the so called self-trapping effect [5]. Above a critical value of the interaction strength, 
the system properties change qualitatively and unbalanced solutions appear, favoring 
one of the wells. A careful discussion of this effect, the relation between mean-field 
and iV-particle behavior as well as its control by external driving fields can be found 
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in [6] . In the many particle system self-trapping is suppressed by tunneling oscillations. 
The same effects are also present in the kicked system. In that case the dynamics of 
the many particle system is much richer. In fact a hole bunch of prominent features of 
quantum dynamical systems like coherent destruction of tunneling and chaos assisted 
tunneling can be found in this model. 

In the following section we present the basic model and give a short review of the 
Floquet formalism. In section [3] we investigate the mean- field dynamics and compare it 
to the much richer many particle dynamics in section HI We conclude the paper with a 
detailed discussion of the tunneling behavior. 



2. Basic model and Floquet formalism 

We focus on a two-site Bose-Hubbard system with the Hamiltonian 

H = | (a\ai - a\a 2 ^j + ^ (a\a 2 + a\a-^j - ^ {a\a\ - a\a 2 ^j , (1) 

where ay, aj are bosonic particle annihilation and creation operators for the jth mode, 
e is the on-site energy difference, v controls the single particle tunneling and c the 
interaction strength. The Hamiltonian commutes with the number operator 

N = a\a\ + a\a 2 (2) 

so that the total number of particles N is conserved. It is convenient to introduce 
angular momentum operators according to the Schwinger representation 

Lx = \ ( a ^ 2 + a ^ ai ) 

Ly = 2~i ( a ^ 2 ~ °^ ai ) ^ 

Lz = \ ( a * ai ~ a ^ a2 ) ' 
which obey the su(2) commutation relations [Li,Lj] = iL^, where k = x,y,z and 
cyclic permutations. The angular momentum quantum number is equal to i = N/2. In 
terms of these operators the Hamiltonian (CQ) assumes the form 

H = eL z + vL x - cL\. (4) 

As already discussed in the introduction, we now consider a periodically kicked 
system, in particular the case of a kicked interaction strength 

+ 0O 

c{t) = ct Sty-mr), (5) 

m=— oo 

i.e. the time-periodic Hamiltonian 

H(t) = H + c(t)V , H = eL z + vL Xl , V = -h\. (6) 

In the symmetric case e = this is in fact the Hamiltonian of a kicked top [14, 15]. A 
similar system with kicked single particle tunneling coupling v has been studied in [17]. 
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Mathematically, the benefit of a kicked system is the enormous simplification of the 
dynamics. The mean-field approximation, for example, can be evaluated in closed form 
as outlined in section [31 

Since the Hamiltonian ([ID is periodic in time, H(t + r) = H(t), it is useful to 
introduce the propagator over one period, the so called Floquet-operator 



F(t) — U(t + T,t) — T 



exp^-i^ H(t')dA 



(7) 



where the operator T takes care of the right time order. The Floquet-operator (JTj) is by 
definition unitary and thus has unimodular eigenvalues 

Ht)> = e-->(t)> . (8) 

The N + 1 quasi-energies e K are the natural counterparts of the energies of time- 
independent systems. We choose e R in the first Brillouin-zone [— 7t/t, 7t/t). 

The eigenstates of the Floquet-operator, the so called Floquet-states \n{t)), form 
an orthonormal basis of the Hilbert space and obey 

|«(t)> = e"^ \<p K (t)) with \<j> K {t + r)> = \<j> K {t)) . (9) 

For the 5-type time dependence the Floquet-operator factorizes 

F = F(0) = ^Llr e -\{eL z +vL x )r_ ^ 

Similar expressions would be obtained by modulating the single particle tunneling 
coupling v and the on-site energy difference e with periodic 5-kicks instead of the particle 
interaction c. In this case one would arrive at a Floquet-operator 

where the order of the two exponential factors of (TlOT) is switched. The stroboscopic 
dynamics induced by iterations of the Floquet-operator, however, is not changed. 

The additional symmetry (in the symmetric case for e — 0) can be expressed via 
an invariance under rotations R x = R x (ii) = e~ mLx around the x-axis about ir. The 
Floquet-operator commutes with R x , 

[F,R x ] = 0, (12) 

and the Floquet states can be separated into two symmetry classes, one even and one 
odd under the rotation R x . If N is even, the angular momentum quantum number I is 
an integer. Thus we have R\ — \ and the eigenvalues of R x are ±1. Equation ^ can 
be reformulated according to 

F\k±) = e- itKT \K±) (13) 

where 

R x \k ± ) = ±\k ± ). (14) 

If iV is odd R x — — 1 holds, since i is half an integer. The eigenvalues of R x then are 
±i. This yields as well two classes of Floquet states \k), which obey 

R x \k±) = ±i |«±) . (15) 

These two symmetry classes will be of importance for the tunneling process investigated 
in the following. 
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3. Mean-field dynamics 

The celebrated mean-field approximation can be achieved by replacing the field operators 
by complex numbers 



a 



-►Vi, 4— ►V'i, J = 1,2. (16) 

Note, however, that the mapping ( |T6l) is not well-defined, since complex numbers 
commute in contrast to operators. Therefore, the operators should be rewritten in 
symmetrized form in a, and a] before the replacement (fT6l) is carried out [10,12]. Here 
the Hamiltonian §J§ is already symmetric, however the number operator ([2]) appears as 

at t t Q l ai + aia l + °2 a 2 + a 2<4 ~~ 2 I / |2 | I / |2 /-, ~\ 

iV = a{ai + a 2 a 2 = — ► +m\ -1 (17) 

and the mean-field wave function is normalized according to l^il 2 + 1^2 1 2 = N + 1 (see 
also [12]). 

This is in fact a classical limit of the system, equivalent to the limit 1/N — > 0, 
where the parameter 1/N takes over the role of an effective Planck constant [15,18]. In 
the case of the Hamiltonian (Tj[|) this leads to a Hamiltonian function 

%,«, 2 ) = f (»1 - V^2) + |(^2 + ^l) - - ^ 2 ^2) 2 (18) 

of the conjugate variables ipj and ip* which obey the classical canonical equations of 
motion 

^ = Wi J = 1 ' 2 - (19) 

The dynamics of such a system can be equivalently put in the form of a discrete 
Gross-Pitaevskii equation (GPE), or discrete nonlinear Schrodinger equation, in terms 
of the nonlinear matrix equation 




e + ck 
v 




(20) 



where K = \ip2 \ 2 — l^il 2 is the population imbalance of the two wells. The dynamics 
described by this equation of motion conserves the norm of the wave function and can 
be conveniently visualized by a mapping onto the Bloch sphere: Similar to ((31), one can 
define the three quantities 

S X = \(^2 + «), S y = 1(^2 - Ml), S z = i(|^l| 2 - |tf ), (21) 

the components of the Bloch vector s G R 3 . The norm of the Bloch vector s — \s\ = 
(N + l)/2 = £+l/2is determined by the normalization of the mean-field wave function. 
Writing the GPE (1201) in terms of the Bloch vector yields the nonlinear Bloch equations 

s £ — £ Sy I 2cs yS % 

Sy = ss x - vs z - 2cs x s z (22) 
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For the periodically kicked system one can, as in the many particle case, define a 
Floquet operator T for the mean-field Bloch vector. Here again the Floquet operator 
factorizes: For c(t) = the dynamics is linear and can be solved in closed form. The 
resulting time evolution can be described by a matrix 

/ ^11 "^r cos cut —-sin cur ^| (1 — cos cur) \ 



K 



cos cur — - smujr 

2 2 

v [1 — cos cut) -sin cut % + \ cos lot 



(23) 



\ ^ I 1— COS LOT ) -SUl CUT ^2+^2 cos CUT / 

a rotation about the axis (v/u, 0, e/uj) with frequency 



u = Ve 2 + v 2 . (24) 

This rotation is followed by a kick, which conserves the ^-component s z , described by 
the torsion matrix 

/ cos2cts 7 sin2cTS 7 \ 



(25) 



— sin 2cts z cos 2cts 2 
V 1J 

The full mean-field dynamics therefore reads 

s m+i = T s m = ICR, s m , (26) 

where s m = s(mr), m G Z. 

In the following numerical studies we will consider a symmetric system (e = 0). 
Stroboscopic iterations of s according to the discrete mapping ( 1261) with v — 1 and four 
different values of the interaction-strength c are presented in figured! Depending on the 
inter particle interaction strength c, the system exhibits the typical behavior of classical 
nonlinear and thus chaotic systems. In the case of no interaction (c = 0) shown in figure 
[TJ a) the system is integrable and one observes the usual Rabi oscillations which give 
rise to simple rotations of the Bloch vector. At the two stable fixed points (±s, 0, 0) the 
population imbalance of the two wells is constant in time. With increasing interaction 
the linear rotation around the x-axis is more and more perturbed by a torsion along 
the z-axis and bifurcations of the fixed points occur. For c = 1.1/(N + 1), shown in 
figure [Hb), the first bifurcation of the fixed point (— s, 0,0) has just caused two new 
fixed points to appear. Increasing c further, the regular islands around the new fixed 
points grow while the fixed points wander toward the poles of the Bloch sphere; the 
example c = 1.8/(A^ + 1) can be seen in figure (He). In addition, chaotic regions appear 
in the vicinity of the unstable fixed point. The two regular islands will be of special 
interest in the following, since they correspond to the famous self-trapping states of 
the non-kicked system as they are localized on the northern and southern hemisphere, 
respectively. Unlike the non-kicked system, self-trapping occurs only in a certain regime, 
since increasing c even further leads inevitably to global chaos. The case c = 3.5/(A^+ 1) 
is shown in figure Q]d); the chaotic sea spreads noticeably over almost the whole sphere, 
sparing only small regular islands. 

The critical value of the parameter c for the occurrence of self-trapping, i.e. the point 
of bifurcation of the prominent fixed point (— s, 0,0), can be calculated in dependence 
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Figure 1. Stroboscopic dynamics of the Bloch vector s for the parameters e — 0, v — 1. 
t = 1 and a) c = 0, b) c = l.l/(N + 1), c) c = 1.8/(JV + 1) and d) c = 3.5/(iV + 1), 
respectively. 



on the other system parameters by linearization. To guarantee stability of the fixed 
point (— s, 0, 0) the interaction strength must obey the inequalities 

=Fl — COSVT ±1 — COSVT . 

— <c< for smvr>0. (27) 

s sin s sin 

Thus, for parameters r — v — 1, this fixed point is stable for c G (—3.661, 1.092)/(A r +l). 



4. Many particle quantum dynamics 

In order to compare the mean-field dynamics with the full quantum dynamics, we 
introduce the SU{2) coherent states, also called atomic coherent states \$,<p). They 
are constructed by an arbitrary 5f/(2)-rotation R($,tp) = exp (i$(L x sin <p — L y cosip)) 
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of the extremal Fock state \N) with all particles in the first well, which is a coherent 
state located at the north pole of the Bloch sphere [18-21], 

\#,<p)=R{#,<P)\N)- (28) 
In the Fock basis they take the form 

= £ \[(V) C0Sn (I) ^ (j)e i(JV - n) * l"> , (29) 

where 

\n) — \n, N — n) — - 1 = (a\f (alf~ n |0, 0) , (30) 
y/n\(J\ — n)\ 

with n = 0, 1, . . . , N, are the usual Fock states. The expectation values of the angular 
momentum operators in a coherent state (Lj) = <p\ Lj <p), j = x, y, z are given by 

N N N 

(L x ) = — sin^cos^, (L y ) = — sin^siny?, (L z ) = — cos-#. (31) 

Thus, the vector (L) = ((L x ) , (L y ) , (L z )) points on a sphere with radius £ = N/2. 

Here again, the Floquet operator induces a discrete mapping from kick to kick. In 
the Heisenberg picture this yields the mapping 

L (m+l) = F f L M F = e i^0T ^icVr L M e - icFr \ e - iHo r ^ 

for the angular momentum operators after the m-th kick. In the symmetric case 
with e = these equations assume the simple form 

F^L X F = ^ [(L x + i(L y cosvr - L z smvT))e- ic(1+2L y sinVT+2LzCOSVT) ] 
+ h.c. 

F^LyF = ^ [(L x + i (L y cos vt - L z sin vr)) e" ic ( 1+2i « sin ^+ 2i - c ° s ^)] 
+ h.c. 

F^L Z F = L y sin-ur + L z cosvr. (33) 

Now we can compare the many particle with the mean-field dynamics. For this 
purpose we focus on the normalized Bloch vector s/s and the normalized vector of 
expectation values (L) /£ of the angular momentum operators. In the classically 
regular regime, the full quantum dynamics closely resembles the mean-field system and 
shows the usual Rabi oscillations according to the single particle tunneling, which are 
modulated by a breakdown and revival scenario due to the discrete spectrum of the 
system. A similar behavior is very common and known for many systems [5,7]. Figure 
[2] shows a comparison of the mean-field and many particle dynamics. Shown is the z- 
component of the angular momentum expectation value (i.e. the population imbalance of 
the two wells) for a symmetric system (e = 0) with v = 1 and a weak interaction strength 
c = 0.5/ (N+ 1), where the initial state is a coherent state \N) located at the north pole, 
or s/s = (0, 0, 1) in the mean-field approximation. Both oscillate with frequency u — 1 
(compare ff2M ) and the envelope of the many particle expectation values decay with a 
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Figure 2. Time evolution of the ^-component of the angular momentum expectation 
value for e — 0, v = 1 and c = 0.5/(iV + 1). The many particle dynamics (• ) for an 
initial coherent state \N) located at the north pole is compared with the mean- field 
trajectory (•) for a starting vector s/s = x = (0,0, 1) (the lines are plotted to guide 
the eye); left: N = 50 particles, right: N = 100 particles. 



width ~ y/N + 1. Note, however, that the deviation between the mean-field and many 
particle dynamics is due to the fact that the mean-field description is based on a single 
trajectory and can be cured to some extent by propagating an ensemble approximating 
initially the coherent state distribution of the many particle system [22]. 

Our interest, however, is not focused on the regular regime, but on the mixed 
regular-chaotic regime with the two prominent regular islands centered at the self- 
trapping states s± = s(i9±,<p±). A mean-field trajectory started in one of these 
states will stay there forever. The full quantum system however, shows a much richer 
behavior here. Starting the propagation in one of the corresponding coherent states 
|=t) = \$±,<P±) the system exhibits dynamical tunneling to the other island, i.e. the 
state |=f). Such dynamical tunneling processes are known for many systems [23,24], 
e.g., a particle in periodically driven double- well potential [25-27] or a periodically 
driven rotor [28-30]. Such a system can be realized with a BEC [31-36] allowing 
direct experimental observation of the tunneling. It should be emphasized, that here 
the dynamical tunneling between both of the island states is a coherent many particle 
process and may not be confused with the single particle tunneling which causes the 
Rabi oscillations. Nevertheless it corresponds to an actual tunneling of the condensate 
between the potential wells, since the northern state is localized almost in the first 
well, whereas the southern state occupies mostly the second well. In figure [3] the 
tunneling process starting from the fixed point on the southern hemisphere in the 
state |^(0)) = |— ) is shown for several parameter sets. An interaction strength of 
c = 2/(N + 1) is chosen to ensure the two regular islands to be not too small and well 
separated. One observes dynamical tunneling through the Bloch sphere. The normalized 
expectation value (L) / i spirals its way on a tube encircling the straight line connecting 
the two island states to the other side of the sphere. The radius of the tunneling tube 
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Figure 3. a) Mean- field dynamics for the parameters £ = 0, v = 1, t = 1 and 
c = 2/(N + 1). the fixed points s± arc marked with red dots, b) Propagation of 
many-particle state |— ) located on the southern fixed point for N = 10 particles; the 
dynamical tunneling to the northern fixed point lasts 40 periods here; c) as before, but 
N = 20 particles; tunneling during 1000 periods d) as in c). but N = 100 particles, 
dynamical tunneling cannot be observed here. 



is determined by the uncertainty of the direction of the angular momentum vector in 
a coherent state, which is proportional to 1/N. The more particles one puts into the 
system, the tighter is the spiral around the connection line of the fixed points and the 
longer gets the tunneling time. However, under some conditions the tunneling time 
can also decrease with a growing number of particles, leading to a strongly enhanced 
population transfer, as will be explained in section |4~T1 But let us first have a closer look 
at the tunneling process. The squared absolute values | (?/;(£) |±)| 2 of the projections of 
\ip{t)} on the coherent island states and the orthogonal subspace for 1^(0)) = |— ) for 
case c) of figure [3] with N = 20 is shown in figure H] as a function of time. After 1000 
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Figure 4. Squared absolute values of the projections of \ip(t)) on the coherent states 
|±) during propagation over 1000 periods for N = 20 particles, e = 0, u = 1, r = 1 
and c = 2/(N + 1); — (— ^>(i)) , — and — orthogonal subspace. 



periods, the population has been transferred almost completely to the state |+) on the 
northern hemisphere. 

For all particle numbers N considered in figure [3j the population |(ac|— )| 2 of the 
southern coherent island state on the Floquet states \k) is significantly large only for 
two of the (N + 1) Floquet states which are denoted by \k±) with corresponding quasi- 
energies e±. As indicated by the ±-sign, they belong to different symmetry classes and 
the island states can be reconstructed up to 94 % by the linear combinations 

|V±> = (|«+> T |«_» (34) 

which hence are a good approximation for the coherent states |±) ~ \ij}±). If one chooses 
as a starting vector, \ip{t)) always stays a linear combination of the \k±) during 
propagation. Remembering equation §§§ for the decomposition of the Floquet states, 
the propagation over n kick periods yields 

|V(nr)) = ^e" ine+T (|« + > + e - in(e — £+)t |k_» (35) 
v2 

= l e - ine+T ((1 - e"^") |+) + (1 + e-' 1 ^) |-)) , (36) 

with (p n = n(e_ — e+)r. The expectation value of the angular momentum varies as 
it i w i-cos^r , l + COS0 n 

{L(nT)} = L + H L_ + sm0 n Im£-_ + (37) 

with L± = (±| L |±), = (— | L \+), where in the present cases ImL h is zero or 

very small so that the ellipse (|3T|) degenerates to a straight line connecting the island 
centers at L_ and L + . Thus the spiraling on a tube around this line is due to the 
other contributing Floquet states. We conclude, that the subsystem consisting of the 
two island states |±) can be described in a good approximation by a two state system 
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Figure 5. Husimi distributions of the island states \ip±); left: Q-{d : <p) = 
|(tf,92|V>-}| 2 , right: Q+(i?,<p) = \(^,f\^+)\ 2 for N = 20 particles, e = 0, u = 1, 
r = 1 and c = 2/(JV + l). 



of the Floquet states |re±). Thus, as any two level system, the system exhibits periodic 
tunneling with period 

2tt 2tt 

-^tunnel = — ~T~ • (38) 

e_ — e + Ae 

These considerations are entirely equivalent to the tunneling process through a 
symmetric potential barrier, where the quasi-energies adopt the function of the energies 
of time independent systems. Thus the tunneling period is solely ruled by the quasi- 
energy splitting Ae. One should keep in mind, that despite the similarities in the 
formalism this is no single particle, but a coherent many particle process. Still the 
dynamical tunneling between the two island states corresponds to an actual tunneling 
of the BEC from one potential well to the other. If the system is localized in the southern 
island |— ), the expectation value of L z is negative and the condensate is located mainly 
in the second well; in the case c = 2/ (TV + 1) over 80% of the condensate is trapped 
there. The opposite is true for the northern island |+). 

To give an impression of the localization of the Floquet states, one can compute 
quantum phase space distributions, e.g. the Husimi distributions [37, 38] 

W,¥>) = l(tf,¥#}| 2 - (39) 
via the ST/ (2) coherent states (p). In figure O the Husimi distribution for the linear 
combinations \ip±) of the two Floquet tunneling states are shown. The sharp localization 
in one of the two regular islands is evident here. 

4-1. Parameter dependencies of the tunneling time 

In general the dynamics of the island states, however, is more complicated and shows 
special features. Tunneling can in fact be extremely stimulated or fully suppressed as 
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Figure 6. Tunneling period T tunnc i calculated from Ae according to ((38)) as a function 
of N for e = 0, v = 1, r = 1 and c = 2/(N + 1); • Ttunnci determined from numerical 
propagation. 

we will see in the following. 

At first we examine the tunneling for different numbers of particles. In figure |6] the 
tunneling period T tunnc i for c = 2/(N + 1) calculated with the help of fl38|) is plotted 
over N; please note the logarithmic scaling of the T-axis. This is in excellent agreement 
with the tunneling time directly determined from the time propagation, plotted in red 
for comparison. For small N the tunneling period increases exponentially. For larger 
N new structures appear, where the tunneling time decreases rapidly about several 
orders of magnitude. These "resonances" have been observed for other systems e.g. the 
anharmonic driven oscillator [25,29,39]. Quite often this phenomenon is called "chaos 
assisted tunneling" [23,25,26,40-42] or "resonance assisted tunneling" [43,44] and can 
be explained as an avoided crossing of the tunneling doublet e± with some third level e c 
which increases the splitting Ae. In such regions one actually has to deal with a three 
state tunneling mechanism. Nevertheless the tunneling period can still be calculated 
via (!38|) quite accurately. 

We illustrate this by means of the resonance at iV « 33 where T tunnel changes by 
two orders of magnitude. The Husimi distributions of the three Floquet states \k±) and 
\k c ) involved in the tunneling for N = 33 are shown in figured The third state \k c ) is 
mostly localized in the chaotic regions of the classical phase space and may therefore be 
denoted as "chaotic" state. The overlap of this state with the island states is actually 
small, |(k c |±)| 2 ~ 0.06. Nevertheless it changes the behavior of the system dramatically. 
To visualize the tunneling process, it is sufficient to focus on the z-component during 
propagation. In figure [8] the ^-component is shown for N = 30 and N = 33. The great 
difference in the tunneling period is obvious here. The system with N = 30 particles 
shows perfect self-trapping up to 10 4 kick periods, whereas the system with just three 
more particles exhibits strong tunneling with T tunne i m 4000r. Thus, the small variation 
of the number of particles obviously changes the entire tunneling process from a two 
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state to a faster three state mechanism. The interplay with the third "chaotic" state 
motivates the name "chaos assisted tunneling" (CAT) for this decrease of the tunneling 
period. If one takes a closer look on the tunneling oscillations of this system the small 
modulations attract attention. This can easily be understood as tunneling into the third 
state |« c ) which is superimposed. The period T c of these modulations can be calculated 
by the means of ( 1581) with the quasi-energy splitting Ae c between the third level and the 
tunneling doublet. In the present case this yields T c = 551 in perfect agreement with 
the numerical propagation. 

To achieve a better understanding of the behavior of the quasi-energies, especially 
crossing or avoided crossing scenarios in regions where the levels really get close to 
each other, one needs to examine their dependence on continuous variables, e.g., the 
kick strength c. Of special interest is the role of the third chaotic state. One has to 
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Figure 8. Population imbalance (L z ) /£ over 10000 kick periods for e = 0, v = 1 and 
c = 2/(N + 1); N = 30 particles, - N = 33 particles. 
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Figure 9. Tunneling time Tt uim ei calculated according to |38|) with the help of Ae 
as a function of c for a fixed number of N = 30 particles, e = 0, v = 1 and r = 1; 
CDT: coherent destruction of tunneling divergences, CAT: chaos assisted tunneling 
resonances, in the regions of the small gaps the two state approximation breaks down. 

keep in mind, however, that the mean-field phase space structure changes by varying c; 
especially both of the regular islands are modified. 

At first we take a look at the tunneling time T tunne i as a function of c where we 
leave the number of particles N = 30 fixed. This plot, shown in figure [9j provides 
similar features as seen for the iV-dependence in figure El For higher values of c there 
exist the same kind of CAT resonances as seen in the N dependence above, where the 
tunneling period decreases by several orders of magnitude. They are due to avoided 
crossing scenarios of the quasi-energy levels, which shall be examined in more detail in 
the following. In figure [10] we show the quasi-energies over the interaction strength c 
for N = 30 particles (note the cylindrical topology; due to the definition of the quasi- 
energies — 7r and 7r have to be identified on the e-axis). Let us focus on the prominent 
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Figure 10. Quasi-energies e K for iV = 30 particles, e = 0, v = 1 and r = 1 as 
a function of c; • tunneling doublet e±; (CDT): real crossing of e± inducing CDT; 
(CAT): avoided crossing with a chaotic level inducing CAT, see text for details. 
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Figure 11. Blow-up of the CAT region of figure [TUl the symmetry of the quasi-energy 
levels are denoted by (±); one observes an avoided crossing of the quasi-energy levels 
with (+)-symmetry; • tunneling doublet e±. 

resonance at c(N + 1) « 2.35 corresponding to the region marked by CAT in the figure. 
A magnification of this region is shown in figure [TT1 Here one of the levels of the doublet 
undergoes an avoided crossing with a chaotic third level e c (since they belong to the same 
symmetry class). The splitting of the avoided crossing is in the same order of magnitude 
as the splitting of the tunneling doublet. Hence, in an intermediate region the three 
states \k±) and \k c ) all take part in a three state tunneling mechanism. The simple two 
state model breaks down here and the tunneling time cannot be evaluated according to 
(158]) (note the small gaps in the graph in figure [9]). Due to the strongly avoided crossing 
the splitting of the tunneling doublet increases, when the third state approaches, which 
causes a significant decrease of the tunneling time. An illustration of this stimulated 
tunneling process, CAT, can be found in [47,48]. 
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Comparing figure M and [TD] it becomes obvious that apart from these resonant 
avoided crossing events there is a large number of avoided crossings which have no 
significant influence on the tunneling period. The avoided crossing between the third 
level and one of the doublet levels has only a very small splitting in these cases, which 
has no influence on the tunneling period, since the quasi-energy doublet stays almost 
undistorted. 

Despite the CAT resonance there is another structure in figure |9] which attracts 
attention: While for small values of c the tunneling period increases faster than 
exponential and shows an almost monotonic behavior, a sharp peak can be observed 
for c pa 1.98/(iV + 1) which in fact is a true divergence. This phenomenon is known 
as "coherent destruction of tunneling" (CDT) [23,45,46]. At these points in parameter 
space the two levels belonging to the tunneling doublet cross. Since they have different 
parity such crossings are not forbidden. This zero splitting automatically leads to a 
divergence in the tunneling period. Thus, the relatively high value of the tunneling 
period that can be observed in figure M for N = 30 particles is due to the direct 
neighborhood to the CDT divergence in parameter space. A diverging tunneling time 
means in fact total suppression of the dynamical many particle tunneling and therefore 
a real self-trapping of the condensate not only in the mean-field approximation, but also 
in the full many particle regime. 

Jf.,2. Tunneling landscapes 

Let us finally explore the whole parameter space available for a symmetric two-mode 
system. One system parameter is the interatomic interaction strength c, the second 
parameter is the time scale of the system, which is given by the single particle tunneling 
coupling v. Both of the effects CDT and CAT are governed by the behavior of the 
quasi-energies in dependence on the system parameters, which in turn is ruled by 
the symmetries of the system. These are mappings under which the eigenvalues and 
symmetry properties of F are preserved; for the Floquet operator ([71) they are given 
by orthogonal matrices. Therefore F has codimension n = 2, which means that at 
least two system parameters have to be varied to enforce a degeneracy of two quasi- 
energies belonging to the same symmetry class [49,50]. Hence accidental degeneracies 
of such quasi-energy planes e K (c,v) occur at most at isolated points in the parameter 
space spanned by c and v. However, planes belonging to different symmetry classes, 
i.e. having opposite parity, actually exhibit exact crossings which form a one dimensional 
manifold in the (c, t>)-plane. Thus, the quasi-energy planes e K (c,v) cross along lines [23, 
sect. 6]. These lines play an important role in the context of tunneling since the 
tunneling period T tunnc i diverges there. Plotting T tunnc i as a function of c and v yields 
the "tunneling landscape" of the kicked double- well BEC. Figure fT2l shows T tunnc i(c, v) in 
false color representation on a logarithmic scale. One observes indeed a rich landscape 
with diversified sequences of ridges and valleys. The sharp lines with small values 
of the tunneling period (dark blue) correspond to the CAT regions; along these lines 
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Figure 12. Tunneling landscape, Xtunnci over (c, v) in false color representation on 
logarithmic scale; N = 30 particles, e = and r = 1. The sharp dark lines with 
reduced tunneling times represent CAT regions whereas these with high values of the 
tunneling time represent CDT regions. 

avoided crossings of the quasi-energy planes with significant splitting occur. On the 
other hand, the sharp lines with high values of T tunne i are actually lines of divergence 
since the quasi-energy planes intersect along these lines in parameter space. Thus, they 
correspond to the CDT regions. In the upper left corner of the presented section of 
parameter space, i.e., for smaller values of c and larger values of v, regular dynamics 
prevails and almost no crossings of the quasi-energy planes take place. The consequence 
of this is a flat tunneling landscape area. The opposite holds in the chaotic regime; here 
a larger number of CAT resonances and CDT divergences occur. Remarkable are the 
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Figure 13. Tunneling landscape, blow-up of figure fT2l avoided crossing of two CDT 
ridges, in the gap one observes a CAT valley. 

"avoided crossings" between both the ridges and the valleys of the landscape. One of 
these is shown in figure fT3l which is a magnification of figure [121 A systematical feature 
at these crossings is the direct neighborhood of resonances and divergences which could 
already be observed in figure [9] that is actually a cut through figure [12] at v — 1. Another 
interesting fact about this tunneling landscape is, that in principle the system can be 
tuned to the CDT regime by adjusting one parameter. Thus for this kicked system 
real self-trapping is possible not only in the mean-field regime, but also for the full 
quantum system and the system can be tuned from enforced tunneling to this regime 
by rather slight variations of the parameters. This yields the possibility of a systematic 
and accurate population transfer between the two potential wells. 

5. Conclusions 

In the present paper we investigated a two-mode BEC with a periodically kicked 
interaction term in a many particle as well as in a mean-field description. This system 
is actually equivalent to the kicked top which is a standard example of quantum chaos. 
Besides the new possible experimental realization of the kicked top, this new context 
sheds light on different aspects of the model. We showed that the tunneling oscillations 
of the many particle system which usually suppress the self-trapping effect show a 
rich spectrum of phenomena like coherent destruction of tunneling and chaos assisted 
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tunneling. The sensitive parameter dependence of the tunneling behavior of the system 
offers an interesting tool for the manipulation of Bose-Einstein condensates. 
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